Quantum Goos-Hanchen effect in graphene 
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The Goos-Hanchen (GH) effect is an interference effect on total internal reflection at an interface, 
resulting in a shift a of the reflected beam along the interface. We show that the GH effect at a 
p-n interface in graphene depends on the pseudospin (sublattice) degree of freedom of the massless 
Dirac fermions, and find a sign change of a at angle of incidence a* — arcsin \j sin a c determined by 
the critical angle a c for total reflection. In an n-doped channel with p-doped boundaries the GH 
effect doubles the degeneracy of the lowest propagating mode, introducing a two-fold degeneracy on 
top of the usual spin and valley degeneracies. This can be observed as a stepwise increase by 8e 2 /h 
of the conductance with increasing channel width. 

PACS numbers: 73.23.Ad, 42.25.Gy, 72.90.+y, 73.50.-h 



Analogies between optics and electronics have inspired 
the research on graphene since the discovery of this mate- 
rial a few years ago [l| . Some of the more unusual analo- 
gies are drawn from the field of optical metamaterials. 
In particular, negative refraction in a photonic crystal 
Q has an analogue in a bipolar junction in graphene if 
the width d of thep-n interface is less than the electron 
wave length \p 3]. Negative refraction is only possi- 
ble for angles of incidence a less than a critical angle 
a c . For a > a c the refracted wave becomes evanescent 
and the incident wave is totally reflected with a shift 
(7 of order Xp along the interface. This wave effect is 
known as the Goos-Hanchen effect Q, after the scien- 
tists who first measured it in 1947. The GH effect was 
already predicted in Newton's time and has become a 
versatile probe of surface properties in optics, acoustics, 
and atomic physics [f|. In particular, the interplay of the 
GH effect and negative refraction plays an important role 
in photonic crystals and other metamaterials @, 0] • 

The electronic analogue of the GH effect has been con- 
sidered previously [j| U [T^, EH , including relativistic cor- 
rections, but not in the ultrarelativistic limit of massless 
electrons relevant for graphene. As we will show here, 
the shift of a beam upon reflection at a p-n interface 
in graphene is strongly dependent on the sublattice (or 
"pseudospin" ) degree of freedom — both in magnitude 
and sign. We calculate the average shift a after multiple 
reflections at opposite p-n interfaces and (contrary to a 
recent expectation [12|) we find that a changes sign at 
a* — arcsin Vsin a c . 

In search for an observable consequence of the GH ef- 
fect we study the conductance of the p-n-p junction, for 
current parallel to the interfaces (see Fig. Q]). (The con- 
ductance for current perpendicular to the interfaces was 
calculated by Pereira et al. [13].) We find that the lowest 
mode in the ra-doped channel has a twofold degeneracy, 
observable as an 8e 2 /h stepwise increase in the conduc- 
tance as a function of channel width. 

We recall some basic facts about the carbon mono- 
layer called graphene [3, EH- Near the corners of the 
Brillouin zone the electron energy depends linearly on 
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FIG. 1: Upper panel: Potential profile of an n-doped chan- 
nel between p-doped regions. Lower panel: Top view of the 
channel in the graphene sheet. The blue solid line follows the 
center of a beam on the A sublattice, while the red dashed 
line follows the center on the B sublattice. The two centers 
have a relative displacement <5o- Upon reflection, each pseu- 
dospin component experiences alternatingly large and small 
shifts a±. 



the momentum, like the energy-momentum relation of a 
photon (but with a velocity v that is 300 times smaller). 
The corresponding wave equation is formally equivalent 
to the Dirac equation for massless spin-1/2 particles in 
two dimensions. The spin degree of freedom is not the 
real electron spin (which is decoupled from the dynam- 
ics) , but a pseudospin variable that labels the two carbon 
atoms (A and B) in the unit cell of a honeycomb lattice. 
To calculate the GH shift we consider a beam, 

* m (z,2/) = J_Jqf{q-qyiy + ^ (^ q( , )/2 J , (1) 

incident on a p-n interface at x — from an n-doped 
region x < 0. The spinor wave function = (\&_|_,\&_) 
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has pseudospin component ty + and ^_ on the A and B 
sublattices. We require that VE ,m is a solution of the Dirac 
equation, 
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U)V = EM, 
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with U = (zero potential in the n-doped region) and 
E = Ep (the Fermi energy). This requirement fixes the 
dependence of the longitudinal wave vector k and the 
angle of incidence a on the transverse wave vector q, 



^(Ep/Tivf 



&rcsm(hvq/Ep). (3) 



For brevity, we will set hv = 1 in some intermediate 
equations (restoring units in the final answers). 

The transverse wave vector profile f(q — q) of the beam 
is peaked at some q € (0, Ep/fiv), corresponding to an 
angle of incidence a — a,rcsm(q/ Ep ) £ (0,7r/2). None of 
our results depend on the shape of the profile, but for 
definiteness we take a Gaussian, 



f(q-q)=cxp{-±(q-q) 2 /A 2 }, 



(4) 



of width A q . 

For Aq small compared to the Fermi wave vector 
kp = Ep/hv we may expand k(q) and u{q) to first or- 
der around q, substitute in Eq. {!]), and evaluate the 
Gaussian integral to obtain the spatial profile of the in- 
cident beam. At the interface x = the two components 
(x exp[-iA^(y - tf£) 2 } of * in (0, y) are Gaussians of 
the same width A y = 1/A g , centered at two different 
mean y-coordinates 



= ±la'{q) = ±\{k F cos a)" 1 . 



(5) 



(The prime in a' indicates the derivative with respect to 
q.) The separation 



S = \y^-y^\ = (k F cos a)' 1 



(6) 



of the two centers is of the order of the Fermi wave length 
Xp = 27r /kp, which is small compared to the width A y 
but of the same order of magnitude as the GH shift - 
so it cannot be ignored. 

Similar considerations are now applied to the reflected 
wave, 



dqf(q - q)e 
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obtained from the incident wave ([T]) by the replacements 
k i— ► —k, a i ► 7T— a and multiplication with the reflection 
amplitude r(q) = |r(<7)|e l ^ 9 ' ) . The two components vj/°_ ut 
of \f out (0,?/) at the interface are Gaussians centered at 

y°± ut = -</>'(<!) T \a\q) = -<f>'{q) =p \{k F cosa)" 1 . (8) 

Comparison with Eq. ^ shows that the first compo- 
nent of the spinor is displaced along the interface by an 



amount cr + = y° ut — y™ = —<fi'(q) — 5q, while the second 
component is displaced by cr_ = y° ut — y™ = —<l>'(q) + So- 
The average displacement, 

<r = \{a + + a-)=-(t>'{q) = -Imj-\nr, (9) 

is the GH shift. As illustrated in Fig. [TJ after N re- 
flections the two components of the spinor are displaced 
by the same amount Na if N is even and by a different 
amount Na =F<5o if N is odd. For TV 3> 1 the difference 2Sq 
between the two displacements becomes small compared 
to the average shift Na. 

The formula ([9]) for the GH shift is generally valid for 
reflection from any interface. To apply it to the step 
function p-n interface we calculate the reflection ampli- 
tude by matching ^ m + <p out at x — to the evanescent 
wave 



dqC(q)e™-^> (*V° {q) . (10) 



K = ^/q 2 -(hv)-*(E F -u y 
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This is a solution of the Dirac equation ^ (with U = Ua 
and E = Ep) that decays into the p-doped region x > 
for hv\q\ > \E F - U \. 

Continuity of the wave function at x = allows us to 
eliminate the unknown function C(q) and to obtain the 
reflection amplitude, 



ie ia (E F -U ) + n + q 



(12) 



E F - U + ie la (n + q)' 

The modulus \r\ = 1 for angles of incidence 

a > a c = arcsin | Uo/Ep — 1| (13) 

such that there is total reflection [16j. Substitution into 
Eq. © then gives the GH shift, 

sin 2 a + 1 - Uo/Ep 
<y = : 

K sm a cos a 
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Xp sin a — sign (Z7 — Ep) sina c 



7r sin 2a 



(14) 
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A negative GH shift (in the backward direction) ap- 
pears at a p-n interface (when Ep < Uq) for angles of 
incidence 



a r < a < a 



arcsin y sin ( 



(15) 



For a > a* the GH shift is positive (in the forward di- 
rection), regardless of the relative magnitude of Ep and 
Uq. In Fig. [2] we have plotted the a dependence of a for 
two representative cases. 

The product a cos a = aj_ is the shift in the direction 
perpendicular to the angle of incidence (while a is mea- 
sured along the interface). This quantity becomes inde- 
pendent of a [in the interval (a c ,7r/2)] when the charge 
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FIG. 2: Dependence on the angle of incidence a of the GH 
shift (j, calculated from Eq. (|14[) for Uo/Ef ~ 1.5 (solid curve, 
p-rt interface) and for Uo/Ef = 0.5 (dashed curve, n-rt in- 
terface). The critical angle for total reflection (below which 
a = 0) equals a c = 30° in both cases. The sign-change angle 
a* = 45° for U /E F = 1.5. 



density in the p-doped region goes to zero at fixed charge 
density in the n-doped region, 



<r± — > l/k F if \E F - U \ < -Bp sin 2 



(16) 



In this limit it does not matter for the sign of the shift if 
Ep is larger or smaller than Uq. Since the perpendicular 
displacement of the two spinor components equals Sq = 
Socosa — 1/kp, the limit (|16p for <jj_ implies that upon 
reflection one component has shift a± — 5^ = equal to 
zero while the other component has shift <j±+Sq = 2/kp. 

As illustrated in Fig.[TJ the GH shift accumulates upon 
multiple reflections in the channel between two p-n inter- 
faces. If the separation W of the two interfaces is large 
compared to the wave length A^, the motion between 
reflections may be treated semiclassically. The time be- 
tween two subsequent reflections is VF/iicosev, so the ef- 
fect of the GH shift on the velocity v\\ along the junction 
is given by 
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wsina + (a/W)v cos a. 



(17) 



Substitution of Eq. (|T4]) shows that, for Uq > Ep, the ve- 
locity vu vanishes at an angle a** satisfying the equation 



• 2 * 

sin a 



(U /E f -1)(kW + 1) 
which for k F W 1 has the solution 
(1 - sina c ) 2 



(18) 



{kp W) 2 sin 2a c sin 2 a c 



0{k F W)- 4 . (19) 



The vanishing velocity shows up as a minimum in the 
dispersion relation, obtained by solving the Dirac equa- 
tion @ with the potential profile 



U(x) 



U for \x\ > W/2, 
for \x\ < W/2. 



(20) 



Matching of propagating waves to decaying waves at x = 
— W and x — produces the following relation between 




FIG. 3: Energy E of waves propagating with wave vector q in 
the y-direction, bounded in the channel — W < x < by the 
potential profile in Fig. [1] The different curves (black solid 
lines) correspond to different modes. (Only the six lowest 
channel modes are shown.) The curves are calculated from 
Eq. (|22[) for UoW/hv = 10 (semiclassical regime). The veloc- 
ity v« — dE/hdq in the y-direction vanishes at the minima of 
the dispersion relation, given by Eq. (|18p (red dashed curve). 
At the (green) dotted lines Tiv\q\ — \E— Uq\ the channel modes 
are joined to modes in the wide region, as indicated schemat- 
ically by the (black) dotted curves. 



E and q: 

[q 2 + E(U - E)] sin kW + kn cos kW 



VE 2 



vV - ((7 - Ef 



(21) 
(22) 



The dispersion relation E(q) is plotted for the first few 
modes in Fig. [3l (A similar dispersion relation was ob- 
tained in Ref. [13J.) The slope determines the velocity, 
Vh = dE/hdq. The minima in the dispersion relation 
where vu — are clearly visible for E ^ Uq. The loca- 
tions of the minima are precisely [TtJ given by Eq. (fl~8"|) 
(red dashed curve). For E > Uq the GH effect increases 
the velocity, which is visible in the dispersion relation as 
a local increase in the slope of the dispersion relation. 
The solid curves in Fig. [3] give the dispersion relation of 
modes that are confined to the narrow n-doped channel. 
At the dotted lines hv\q\ = \E — Uq\ these channel modes 
are joined to the modes in the wide p-doped region (as 
indicated by the dotted curves in Fig. [3]). 

As the channel width is reduced so that UqW/Tlv be- 
comes of order unity, we enter the fully quantum me- 
chanical regime. The minimum in the dispersion relation 
becomes very pronounced for the lowest channel mode, 
as we show in Fig. |4l There are two minima at q « 1/W 
and q as —1/W, each contributing to the conductance a 
quantum of e 2 /h per spin and valley degree of freedom. 
The total contribution to the conductance from the low- 
est channel mode is therefore 8e 2 /h. As shown in Fig. 
[SJ if W is reduced further, the two degenerate minima in 
the dispersion relation merge into a single minimum at 
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FIG. 4: Same as Fig. [3] but now showing the lowest channel 
modes in the fully quantum mechanical regime UoW/hv — 3. 
The two minima at q — ±0.83 W~ l each contribute indepen- 
dently an amount of 4e 2 //i to the conductance. 
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FIG. 6: Conductance versus channel width, calculated nu- 
merically at zero temperature (thin red curve) and at a fi- 
nite temperature (thick black curve) . The dashed black curve 
gives the number of propagating modes, calculated from the 
dispersion relation. 
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FIG. 5: Plot of the lowest mode for three values of UoW/hv, 
showing how the the two minima merge into a single minimum 
at q = upon reducing W . 



q = (this happens at UoW/hv — 1.57), and for smaller 
W the lowest channel mode again contributes the usual 
amount of 4e /h to the conductance. 

To test these analytical predictions, we have performed 
numerical simulations of electrical conduction in a tight- 
binding model of a graphene sheet covered by a split- 
gate electrode. The geometry is similar to that studied 
in Ref. [l8j (but not in the p-n junction regime of interest 
here). Using the recursive Green function technique on 
a honeycomb lattice of carbon atoms (lattice constant a) 
we obtain the transmission matrix i, and from there the 
conductance G = (2e 2 / 7i)Tr trf . Only the twofold spin 
degeneracy is included by hand as a prefactor, all other 



degeneracies follow from the simulation. The graphene 
strip is terminated in the re-direction by zigzag bound- 
aries (separated by a distance Wtotai = 220 a), while it 
is infinitely long in the y-direction. A smooth poten- 
tial profile defines a long and narrow channel of length 
L = 1760 a and a width W which we vary between 
and 30 a. The potential rises from in the wide reser- 
voirs (far from the narrow channel), to Uq — 0.577 Tiv/a 
underneath the gate, and has an intermediate value of 
^channel = 0.277 hv / a inside the channel (where the gate 
is split). The Fermi energy is kept at Ef — 0.547 Tiv/a, 
so that it lies in the valence band underneath the gate, 
while it lies in the conduction band inside reservoirs and 
channel. 

Results of the simulations are shown in Fig. [BJ From 
the dispersion relation we read off the total number of 
propagating modes (dashed curve). The zigzag edges 
of the graphene strip support one spin-degenerate edge 
mode, so the conductance levels off at 2e 2 /h as the chan- 
nel is pinched off. Upon widening the channel, the new 
channel modes have the 8-fold degeneracy predicted by 
our analytical theory. The valley degeneracy is not exact 
(notice the small intermediate step at W = 20 a) , as ex- 
pected for a finite lattice constant. The zero-temperature 
conductance (thin red curve) shows pronounced Fabry- 
Perot type oscillations, due to multiple reflections at the 
entrance and exit of the channel, with an envelope that 
follows closely the number of propagating modes. 

At finite temperature (black curve) the oscillations are 
averaged out, but the excess conductance characteristic 
of the Goos-Hanchen effect remains clearly observable 
at the temperature T = 0.02(?7o — Ep)/kB used in the 
simulation. Scaling up to realistic parameter values, we 
can set the channel width W = 100 nm at the first con- 
ductance step, hence U — E F = 0.03hv/W ~ 10 K, so 
this would correspond to a temperature of 0.2 K. The 
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Fermi wave length Xp in the channel is of order 100 nm 
for these parameter values (of the same order as W at 
the first step) , well above the typical width d ~ 40 nm of 
a, p-n interface [19| . Note that d is two orders of magni- 
tude larger than a = 0.25 nm, so the potential is indeed 
smooth on the scale of the lattice constant (as assumed 
both in the analytical and numerical calculations). For 
ballistic transport through the constriction the mean free 
path should be well above the 100 nm scale. 

In conclusion, we have identified and analyzed a novel 



pseudospin-dependent scattering effect in graphene, that 
manifests itself as an 8e 2 /h conductance step in a bipolar 
junction. This quantum Goos-Hanchen effect mimics the 
effects of a pseudospin degeneracy, by producing a pro- 
nounced double minimum in the dispersion relation of an 
n-doped channel with p-doped boundaries. Such a chan- 
nel can be created electrostatically, and might therefore 
be a versatile building block in an electronic circuit. 
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Foundation NWO/FOM. 
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